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[l] We exploit quantitative metrics to investigate the 
information content in retrievals of atmospheric aerosol 
parameters (with a focus on single-scattering albedo), 
contained in multi-angle and multi-spectral measurements 
with sufficient dynamical range in the sunglint region. The 
simulations are perfonned for two classes of maritime 
aerosols with optical and microphysical properties compiled 
from measurements of the Aerosol Robotic Network. The 
information content is assessed using the inverse formalism 
and is compared to that deriving from observations not 
affected by sunglint. We find that there indeed is additional 
information in measurements containing sunglint, not just 
for single-scattering albedo, but also for aerosol optical 
thickness and the complex refractive index of the fine 
aerosol size mode, although the amount of additional 
information varies with aerosol type. Citation: Ottaviani, M., 
K. Knobelspiesse, B. Cairns, and M. Mishchenko (2013), 
Infonnation content of aerosol retrievals in the sunglint region, 
Geophvs. Res. Lett., 40, 631-634, doi:10.1002/grl.50148. 

1. Introduction 

[ 2 ] While the radiative forcings exerted by some atmo- 
spheric constituents (such as greenhouse gases) are known 
with satisfactory accuracy, the role of aerosols bears uncer- 
tainties so large that they prevent climate models from 
running at the desirable accuracy [Hansen et al., 2011; 
Penner et al., 2011 \Loeb and Su, 2010]. The impact of aero- 
sol particulates is indeed complex and includes both direct 
and indirect effects, leading to positive and negative forcings 
via a myriad of feedback mechanisms. Direct effects are 
mostly determined by the microphysical properties, which 
govern the ratio between absorption (leading to positive 
feedback or wanning) and scattering (leading to cooling 
when radiation is reflected back to space) of electromagnetic 
radiation. Highly reflecting sulfate aerosols can originate 
from volcanic eruptions [Bates et al., 1992] and fossil fuel 
combustion [Smith et al., 2001] or within the silicates which 
predominantly contribute to the composition of dust raised 
by sandstonns [Wagner et al., 2012], Carbonaceous particu- 
lates, originating from biomass, biofuel, and fossil fuel burn- 
ing, are characterized instead by non-negligible absorption 
[Moosmiiller et al., 2009], and much emphasis has recently 
been placed on the uncertainty associated with the efficiency 
of this specific process [Andreae, 2001], 
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[ 3 ] The parameters of importance in aerosol retrieval are 
the column optical thickness, the effective radius and vari- 
ance, and the complex refractive index. The typical bimodal 
nature of aerosol populations requires these parameters to be 
determined for both modes. The overall situation is further 
complicated by the extensive variability of aerosol distribu- 
tion, deriving from regional emission sources and from the 
vertical assortment linked to the dynamics of transport pro- 
cesses [Seinfeld et al., 1998]. 

[ 4 ] A decade ago, Kaufman et al. [2002] hypothesized that 
sunglint, the strong signal caused by reflection of sunlight 
from water surfaces, could be exploited to improve the retrie- 
vals of aerosol absorption. The strategy envisioned the exploi- 
tation of off-glint regions to constrain the scattering properties 
of the aerosol, especially feasible if polarimetric measurements 
are available [Mishchenko and Travis, 1997], together with 
direct transmittance measurements at the center of the glint, 
where the higher signal-to-noise ratio would arguably benefit 
the estimate of extinction. Despite the efforts to improve the 
description of the sunglint phenomenon [Kay et al., 2009], 
and include it in radiative transfer codes as a boundary condi- 
tion [Ottaviani et al., 2008], a rigorous assessment of the fea- 
sibility of improved retrievals in glint regions is still missing. 
Because of this, and the lack of sufficient dynamical range, a 
wealth of satellite observations is systematically discarded. A 
clear example concerns the ocean color community, since each 
composite image of global chlorophyll concentration exhibits 
periodic, wide latitudinal swaths of missing data. 

[ 5 ] Inverse methods are most general and can be applied to 
the widest class of retrievals. Well-established inversion 
schemes [Rodgers, 2000], to predict the uncertainty of param- 
eters retrieved from an observation of given characteristics 
(spectral range, radiometric uncertainty, etc.), have been used 
to test the retrieval capability of multi-angle, multi- and hyper- 
spectral, and polarimetric instruments [Rodgers and Connor, 
2003; Lebsock et al., 2007; Hasekamp and Landgraf 2007; 
Waquet et al., 2009; Hasekamp, 2010; Knobelspiesse et al., 
2011; Coddington etal., 2012]. Adhering to the implementation 
of Knobelspiesse etal. [2012], who compared contemporary in- 
strument designs for the retrieval of fine mode aerosol proper- 
ties, we investigate the potential value of measurements taken 
within the sunglint region. 


2. Methodology 

[6] The act of measurement always introduces some 
noise e in the relation linking the //-dimensional state 
vector of a system x to an rw-dimensional measurement 
vector y = F(x) + e, where F is the forward model which 
describes the knowledge of the measurement process 
and the physics of the problem (which may be imper- 
fect). The forward model can often be locally linearized 
about a reference state x 0 : 
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y-F(xo) =^^(x-x 0 ) + e = K(x-x 0 )+e , (1) 

with the m x n Jacobian matrix K i j=dF i (x)/dxj expressing 
the sensitivity of the model to an infinitesimal change in 
each parameter. Each of the m rows of K corresponds (i.e., 
maps) to a coordinate in measurement space. 

[ 7 ] To know how a measurement updates a certain prior 
knowledge P{x) of a state, we can use Bayes’ theorem which 
provides information on how the measurement probability 
density function (pdf) maps into state space and combines 
with the a-priori pdf to form an a-posteriori pdf / > (x|y) for 
the state: P(x\y) = P(y\x)P(x)/P(y). The probability P(y|x), 
of obtaining y if the state is x, is described by the statistics 
of the measurement noise since x is mapped into the region 
of measurement space enclosed by the measurement error. 
P( y) is often disregarded as a normalization factor. For linear 
problems with Gaussian statistics, Bayes’ theorem enables 
the derivation of expressions for the expected value x of 
the state vector and for the associated n x n, retrieval error 

A 

covariance matrix S: 

x = x„ + (K T S-'K+ S“ I ) _I K r S“ 1 (y - Kx a ) , (2) 

=K r S ( T 1 K+S fl 1 (3) 

with the superscript 1” indicating the matrix inverse and “T” 
the matrix transpose. Equation (3) expresses the connection 
between instrument characteristics and retrieval uncertainty. 
The m x m, symmetric and non-negative measurement error 
covariance matrix S e describes the uncertainties associated 
with each measurement and their correlations. The n x n, 
a-priori error covariance matrix S„ measures the uncer- 
tainties in x a , a guess of the state prior to making the 
measurement. Ultimately, the minimum of the cost func- 
tion i 2 = (y - Kx) r S ( T 1 (y - Kx) + (x - x a ) r S“ 1 (x - x a ) 
is found at x, where the expected value of / 2 should be equal 
to the number of degrees of freedom m if the model is valid. 
When a-priori information is absent, the expression reduces 
to that for a weighted least squares regression. 

[s] The Shannon information content (SIC), 

H = — (log 2 |S _1 S fl | where | • | is the determinant, is readily 
available once the covariance matrices are computed. This 
quantity expresses the gain of knowledge in bits of informa- 
tion consequent to the act of measurement, and it bears strict 
analogy with the thermodynamic definition of entropy. In state 
space, it measures the reduction of the volume enclosing the 
accessible states due to the observation. 

[ 9 ] The simulations in this study were devised for a hypo- 
thetical instrument that can be considered as a common denom- 
inator among those already available in orbit. We examined 
along-track measurements of total reflectance only, at seven 
viewing angles as available with the Multi-angle Imaging 
Spectral Radiometer (MISR; [Diner et al., 1998]) and the 
Polarization and Directionality of the Earth’s Reflectances 
(POLDER; [Tanre etal., 2011]). The MISR sensor has actually 
nine cameras, but observations at ±70° were excluded because 
of concerns that vertical structure could create ambiguities at 
such shallow geometries. We also limited the number of avail- 
able wavelengths to five channels in the visible (410, 440, 560, 
670, and 870 mn) as POLDER and MISR do not have channels 
in the Short-Wave Infrared (SWIR). 


[ 10 ] The resulting measurement vector contains a set of ra- 
diance measurements Rj made at various viewing angles (9,- 
for each spectral channel f, y= {RfO i,Ai), ■ ■ .,Rj(6 m ,l m )}. 
It is advantageous to use dimensionless reflectance units to 
keep the signal from different channels on the same scale. 
The relationship to intensity / is Ri = (/7I7q) / (FoCOsGj) , 
where r 0 is the solar distance in astronomical units used to 
obtain the proper irradiance from the annual average exo- 
atmospheric irradiance F 0 (W m~ 2 ), and Q s is the solar zenith 
angle. The measurement uncertainty associated with each Rj 
is composed of the instrument dark current noise ay-, the 
shot noise c> a , and the calibration uncertainty a h \ 

oi, = 2 ofiRi/I) 2 + afRj/I) + ciR) . (4) 

[ 11 ] In this work, we assume a standardized model of in- 
strumental uncertainty [Knobelspiesse et al., 2012] with 
0y=7 x 10“ 5 , (J a = l x 10“ 8 , and ff A = 0.03. 

[ 12 ] Equation (4) defines the diagonal elements 
S C jj = ( a 2 Ri )i (for i = j). Uncertainties due to dark current 
and shot noise are assumed to be random, and thus uncorre- 
lated. A description of angle-to-angle correlation in calibra- 
tion errors, normally neglected for simplicity, is attempted 
as in Knobelspiesse et al. [2012] with a Markov process 
generating off-diagonal elements with the form: 

S t .ij = i fj, same channel, _ 

S Ci ij = 0 if j , different channels, 

where a is the angular separation between adjacent, instanta- 
neous fields of view in one scan, and p is a correlation param- 
eter ranging from p = 0 (totally uncorrelated system) to p = 1 
(fully correlated system). Here, we assumed p = 0.9 to repre- 
sent the high angular correlation typical of calibration errors. 

[ 13 ] The forward engine at the core of the procedure is based 
on the doubling-adding method [De Haan et al., 1 987]. The lin- 
earization process is repeated for eight different aerosol optical 
thicknesses (AOTs) between 0.039 and 1.0 at 560 nm in order 
to cover a large portion of the space of plausible aerosol loads. 
For each AOT, the Jacobians are estimated as finite derivatives 
via perturbing each state parameter xj by a small amount. From 
the widely used suite compiled by Dubovik et al. [2002], and 
based on the large statistical sample of AErosol RObotic 
NETwork (AERONET; [Holben et al., 1998]) observations, 
we selected two aerosol classes representative of a maritime 
scenario (see Table 1): the non-absorbing Lanai (LH) type, 
for which the AOT is equally partitioned between fine and 
coarse modes, and the moderately absorbing Maldives (MI) 
type, for which the fine mode is responsible for 90% of the 
AOT and the coarse mode for the remaining 10%. The aerosols 
are assumed to consist of homogenous spheres and are placed 
in the first kilometer of a parallel-slab atmosphere to mimic a 


Table 1 . Aerosol Microphysical Properties Used in This Study 3 



Maldives 

Lanai, HI 

Refractive index 

1.44-0.01 10i 

1.36-0.0015/ 

reff,f 0<m), v e ff f 

0.222, 0.236 

0.201, 0.259 

reff.c (/an), v efEc 

4.96, 0.782 

4.29, 0.588 

AOTf(560 nm)/AOT tot (560 nm) 

0.9 

0.5 


“From the AERONET climatology in Dubovik et al. [2002], The spheri- 
cal, homogenous particles have a spectrally invariant refractive index. Size 
is specified bimodally with lognonnal distributions. 
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standard boundary layer. The climatology of all aerosol classes 
is the a-priori knowledge used to compile the matrix S a . 

[ 14 ] As with virtually every radiative transfer code capable 
of including a rough ocean as a lower boundary condition, 
we employ the model by Cox and Munk [1956], which para- 
meterizes the wave slope distribution as a linear function of 
wind speed. The simulations were run for a solar zenith an- 
gle of 45° and a low wind speed (3 m s -1 ) so that the glint is 
bright and limited in extent. This allows us to consider, as a 
first scene, scanning observations taking place along the 
principal plane of reflection where the sunglint signal 
reaches its peak. The second scene can then be selected at 
a relative azimuth of 30°, where sunglint is nearly absent 
and the range of available viewing angles (which for scan- 
ners depends on the orientation relative to the principal 
plane) is only moderately affected. The underwater portion 
of light scattering was modeled with a value of chlorophyll 
a concentration of O.lmgm -1 , representative of standard 
open-ocean waters [Chowdhary et al., 2012], 

3. Results and Discussion 

[ 15 ] The absolute uncertainties as a function of AOT associ- 
ated with retrievals from in-glint observations are depicted 
in Figure 1 with orange lines, and those from off-glint observa- 
tions with blue lines. The solid and dashed line styles are for 
the Ml and LH aerosol types, respectively. The values in the 
top-right comer of each panel represent the starting (a-priori) 


value for the specific parameter: uncertainties remaining close 
to this initial choice are indicative of poor retrieval capabilities 
for the associated quantity. This kind of behavior is shown, for 
example, by the effective variance of the fine mode, which is 
notoriously difficult to determine especially if measurements 
of polarization are not available [Knobelspiesse et al., 2012], 

[ 16 ] Increased accuracy in the determination of the total 
optical thickness from glint observations stems only in part from 
a better estimate of absorption. Regarding the long-standing 
question concerning the Single-Scattering Albedo (SSA), the 
detectable improvement (-0.01) observed for moderate AOTs 
of the absorbing MI aerosol (solid lines) degrades at high 
AOTs. The larger improvement for the LH aerosol is deceptive, 
since this aerosol type is minimally absorbing. The same argu- 
ment is valid for the imaginary part of the refractive index of 
the fine mode from which the SSA is derived. 

[ 17 ] The relative magnitude of the uncertainties for the 
parameters linked to the coarse mode of the two aerosol 
models is not surprising, because for the MI type the fraction 
of this mode is much smaller. When the coarse mode load is 
significant, these uncertainties can be greatly reduced by 
including channels in the SWIR. 

[is] An accurate retrieval of the wind speed clearly depends 
on the availability of a glint profile to fit the Cox-Munk model 
(i.e., it benefits from scans oriented closely to the principal 
plane) and degrades at higher AODs because the sunglint 
patch becomes partly obscured. The same reasoning applies 
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Figure 1 . Simulated absolute uncertainties (except for the bottom-right panel, representing the SIC) as a function of optical 
thickness for the retrieval of relevant parameters of non-absorbing Lanai (LH, dashed lines) and absorbing Maldives (MI, solid 
lines) maritime aerosols (see Table 1). Results from observations within the glint are in orange, and those from off-glint scans in 
blue. At the top-right comer of each panel is the a priori value for the relative parameter. The asterisk indicates a derived param- 
eter, not directly retrieved, and for which the uncertainty is calculated using the uncertainties in its constituent parts. 
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to the determination of the chlorophyll concentration from the 
isolation of the water leaving radiance contribution, i.e., the at- 
mospheric correction is more accurate at low optical thick- 
nesses. Moreover, in this case the in-glint and off-glint results 
nearly overlap since the glint is an interface effect that has no 
influence on the underwater field (provided that sufficient dy- 
namical range is available). Although the sensitivity to the ver- 
tical distribution of aerosols is also improved, the uncertainty 
remains high (~ 50%), which confirms the challenges encoun- 
tered by passive sensors in determining such a parameter, es- 
pecially in the presence of absorption. 

[ 19 ] A glance at the SIC suggests that glint measurements 
of moderate-to-high aerosol loads add 6 to 10 bits of infor- 
mation compared to off-glint observations, depending on 
the aerosol type. For the absorbing MI model, the measure- 
ment provides the capability of discerning up to 2 6 ~ 60 ad- 
ditional states relative to the respective a priori number. 

[20] The formalism employed here fully predicts the ac- 
tual performance of a sensor only if all sources of measure- 
ment and model uncertainties are properly accounted for. 
The results quantify the sensitivity to parameter variability, 
but should not be taken as an absolute measure of parameter 
retrievability. Considering that the physics of the system 
cannot be perfectly captured by the forward model, and that 
sub-pixel variability is neglected, the uncertainties should be 
taken as underestimates. Moreover, the ideal form of the 
measurement error covariance matrix greatly depends on 
instrument-specific correlations. Nevertheless, the objec- 
tive of this study is not affected by these limitations, since 
in- and off-glint observations differ exclusively in the 
strength of the illumination reflected by the surface. 

4. Conclusions 

[21] This work focused on sensitivity of simulated, remote 
sensing observations of tropospheric aerosols to absorption. 
The analysis pertains to maritime aerosol models compiled 
after AERONET measurements and supports the hypothesis 
of Kaufman et al. [2002] that the increased signal-to-noise 
ratio, obtained by scanning across the glint region, can in- 
deed improve the retrieval of aerosol absorption. The effect 
is nevertheless moderate and is a consequence of the reduc- 
tion in uncertainty of all descriptive parameters. 

[22] We have chosen to model the prototype instrument for 
which the simulations are performed with characteristics sim- 
ilar to those of sensors already in orbit. It should be noted that 
the computed uncertainties undergo expected and very signif- 
icant reductions by including channels in the SWIR and polar- 
ization capabilities. In the first case, improvements are noted 
mostly in the uncertainties relative to the aerosol coarse mode. 
In the case of polarization, it is clear that the exquisite sensitiv- 
ity to particle microphysics benefits the retrieval of both fine 
and coarse aerosol parameters, to the extent that the presence 
or absence of the glint becomes nearly irrelevant. 

[23] Acknowledgments. This work is dedicated to the memory of Yoram 
Kaufman. Partial support from the Gloiy Mission Project and the Radiation 
Sciences Program managed by Hal Maring is gratefully acknowledged. 
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